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ABSTRACT 

We report the discovery of a supernova (SN) with the highest apparent energy output to date 
and conclude that it represents an extreme example of the Type Iln subclass. The SN, which was 
discovered behind the Large Magellanic Cloud at z = 0.289 by the SuperMACHO microlensing survey, 
peaked at Mr = —21.5 mag and only declined by 2.9 mag over 4.7 years after the peak. Over this 
period, SN 2003ma had an integrated bolometric luminosity of 4 x 10 51 ergs, more than any other SN 
to date. The radiated energy is close to the limit allowed by conventional core-collapse explosions. 
Optical spectra reveal that SN 2003ma has persistent single-peaked intermediate-width hydrogen 
lines, a signature of interaction between the SN and a dense circumstcllar medium. The light curves 
show further evidence for circumstcllar interaction, including a long plateau with a shape very similar 
to the classic SN Iln 1988Z - however, SN 2003ma is ten times more luminous at all epochs. The 
fast velocity measured for the intermediate-width Ha component (~6000 kms -1 ) points towards an 
extremely energetic explosion (> 10 52 ergs), which imparts a faster blast-wave speed to the post-shock 
material and a higher luminosity from the interaction than is observed in typical SNc Iln. Mid-infrared 
observations of SN 2003ma suggest an infrared light echo is produced by normal interstellar dust at 
a distance ~0.5 pc from the SN. 

Subject headings: circumstellar matter — supernovae: individual (SN 2003ma) — dust, extinction 
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1. INTRODUCTION 

Hydrogen-rich (Type II) core-collapse supcrnovae 
(SNe) from massive stars typically have peak luminosi- 
ties corresponding to absolute magnitudes — 14 > Mr > 
— 18 mag that are powered by thermal energy that is 
deposited into the expan ding SN envelope during shock 
breakout (|Li et al.H2010t ). Depending on the mass of the 
hydrogen (H) envelope, a plateau can arise in the op- 
tical light curve from a cooling wave of H recombina- 
tion that recedes through the ejecta layers. After this 
photospheric phase, the light curve displays an exponen- 
tial decline when it is predominantly powered by the ra- 
dioactive decay of isotopes created in the explosion ( 56 Ni 

In the presence of a dense circumstellar medium 
(CSM), the expanding SN shock will collide with the 
CSM and con vert the bulk kinetic energy of the explo- 
sion into light (jChevalier fc Fransson 1994J), and produce 
a relatively narrow (2-4 xlO 3 kms" 1 ) Ha line that 
is the hallmark of SNe Iln (ISchlegell [199(1 . The addi- 
tional luminosity from the circumstellar interaction can 
be extremely large and even dominate the total luminos- 
ity. Furthermore, while the luminosity from radioactive 
decay will decline quickly after a few months, the lu- 
minosity from circumstellar interaction can persist at a 
near-constant level for several years. 

Recently, several extremely luminous SNe have 
emerged with peak luminosities with Mr < —21 mag. 
SN 2006gy had a very lo ng rise time (70 days) and peaked 
at M r = -21.8 mag ([Smith et all I2007U lOfek et al.l 
120071) . SN 2008fz had a similar pea k brightness 
(My = —22.3 mag) and light-curve shape ([Drake et al.l 
[20T1 . SN 2005ap was very luminous at peak (M un f = 
—22.7 mag), but had a fast rise and decay. SN 2008es 
was similar to SN 2005ap, peaking at My = —22.3 ma ; 
and having a fast rise and decline (jMiller et al.l |200~ 
iGezari et al.l l2009a|) . A self-ob scured luminous super - 
nova, SN 2007va, was detected bv lKozlowski et al.l (|2010l ) 
only in the infrared with an absolute mid -IR peak mag- 
nitude of M[ 4 . 5 ] —24.2. More recently, iQuimbv et al.1 
(|2009f ) and lGal-Yam et al.l (120091 ) announced the discov- 
ery of three extremely luminous SNe, respectively, which 
the y suggest are po wered by a pulsational p air instabil- 
ity (jWooslev et all 120071) . iPastorello et al.l (|2010l ) finds 
evidence that these extremely luminous SNe may be con- 
nected to SN Ic. 

A pulsational pair instability is expected to be impor- 
tant in only the most massive stars - those exceeding 
95 Mq. The production of electron-positron pairs re- 
sults in a contraction and then explosive nuclear burn- 
ing which ejects some significant number of solar masses 
worth of material from the envelope. Subsequent repe- 
titions of this sequence of events result in ejected shells 
catching up with previous ejected material, now at much 
larger radius, producing radiated energy due to the shell 
collisions. This process is estimated to be capable of pro- 
ducing 10 50 ergs of light and, as importantly, can repeat 
on short timescales, providing longer-timescale luminos- 
ity. 

The energetics of these events push the envelope of 
our understanding of stellar evolution. The peak lu- 
minosity, if powered by radioactive decay, would re- 
quire ~10 Mq of 56 Ni. This very large amount points 



to an extremely massive progenitor and the possibility 
that these events w e re the result of a pair-instability SN 
(|Barkat et all 119671 ; iBond et al.l I1984T) . Alternatively, a 
significant amount of the luminosity may be produced 
by circumstellar interaction, but the mass of the CSM 
necessary for the luminosity still indicates that the pro- 
genitors must have been massive stars with significant 
mass-loss histories. 

Furthermore, there is a limit on the energy for which 
conventional core-collapse explosions are viable. Under 
a conventional core-collapse explosion, the maximum en- 
ergy emitted by a SN is equivalent to the rest mass of a 
neutron star, a few 10 54 ergs, with 99% being emitted as 
neutrinos (e.g.. iWooslev fc Ja nka 2005|). The remaining 
energy, a few 10 52 ergs, is either coupled to the baryonic 
material as kinetic energy, or emitted as electromagnetic 
radiation. If a SN has demonstrably greater than a few 
10 52 ergs, the conventional core-collapse scenario must 
be rc-cxamined. This argument has been used for the 
extremely energetic broa d-lined SNe Ic associated with 
gamm a-ray bursts (e.g., iWoosleyl 119931 : llwamoto et al.1 
119981 ) ; although these constraints are placed on this class 
of objects using the kinetic energy, not the radiated en- 
ergy. 

In the local universe, there ha ve been several wel l- 
observe d SNe Iln: SNe 1988 Z dTuratto et all fl993T) 
1994 W (|Sollerman e Fall fl998h . 1995N dFransson et al.l 
I20021). 199 8S (e.g.. iLeonard et al.l 120001). 1 999el 
( Pi Carlo et all 120021) 2005ip (|Smith et al.ll2009bD . and 
2007rt (|Trundle et al.1120091 ). In addition to their simi- 



lar spectral evolution, SNe Iln all have long plateaus in 
their light curves after an initial decline. This plateau 
is from shock energy being continuously converted into 
visual light through circumstellar interaction. Once the 
shock extends to a radius where the density of the CSM 
drops, the source of the luminosity is diminished and the 
SN fades significantly. 

Observations suggest that dust can form in the 
ejecta of SNe dMoselev et al.lll989UKozasa et all 119891: 



Sugerman et al. 



20061: iSmith et al.M2008bl : iRho et atl 



20081: IFox etHf 120091: iKotak et alll2009D . This process 



appears to be enhanced by increasing the density of the 
CSM; i.e., SNe Iln appear to be better at forming dust 
than SNe IIP. An infrared light echo can appear at late 
times when light from the SN explosion reaches a dust 
shell, and th e UV/optical light is reprocessed into the 
infrared fe.g.. iGerardv et al"1l2002aD . 

Here we present SN 2003ma, discovered behind 
the LMC by the SuperMACHO microlcnsing survey. 
SN 2003ma has the light curve shape and spectral sig- 
natures of a SN Iln, but with a peak luminosity compa- 
rable to the most luminous SNe ever discovered. In Sec- 
tion [3 we describe the densely-sampled SuperMACHO 
and OGLE-III difference imaging photometry (covering 
a baseline of 12 years, and the SN light curve over 4.7 
years) , and spectroscopy taken at the peak of the SN, and 
1, 2, 5, and 6 years later, and in Section [3] we compare 
the observations to well studied SNe Iln. In Section @] 
we calculate the total bolometric luminosity and energy 
output of the SN, and show that it exceeds the radiative 
output of any other SN observed to date. We compare 
possible models for the source of the energy, and show 
evidence for an IR echo from variable IR flux measured 
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Fig. 1. — SuperMACHO images of the region around SN 2003ma. 
Left: Template image from 2001 November 20 (pre-event). Mid- 
dle: Search image from 2003 December 17, showing a detection of 
SN 2003ma at the center of the image. Right: Difference image 
(center minus left with an inverted color scale) clearly showing the 
detection of SN 2003ma. The white line in the lower right corner 
indicates 2". North is up and East is left. 



in archival Spitzer IRAC observations. In Section [5l we 
summarize our results. 

2. OBSERVATIONS & REDUCTIONS 
2.1. Optical Photometry 

Starting in 2001, the SuperMACHO Project microlens- 
ing survey used the CTIO 4 m Blanco telescope with 
its 8K x 8K MOSAIC imager (plus its atmospheric dis- 
persion corrector) to monitor the central portion of the 
LMC every other night for 5 seasons (September through 
December). The images were taken through our cus- 
tom "VR" filter (A c = 625 nm, SX = 220 nm; NOAO 
Code c6027) with exposure times between 60 and 200 
seconds, depending on the stellar density of each field. 
In addition to the VR filter, images were occasionally 
obtained through the B and I filters. Throughout this 
paper, we denote brightnesses measured in the natural 
CTIO ma gnitude system a s B^m , VRsm and Iqm as 
defined in IMiknaitis et all (|2007T ) . To search for vari- 
ability, PSF-matche d template images were subtracted 
from search images jRestet all 120051: iGarg et all 120071: 
IMiknaitis et al.|[2007D . The resulting difference images 
are remarkably clean of the (constant) stellar background 
and are ideal for searching for variable objects. Our 
pipeline detects and catalogs all variable sources. 

On 2003 December 13, the SuperMACHO survey de- 
tected a non-microlensing transient event at position 
RA = 05:31:01.878, Dec = -70:04:15.89 with a VR mag- 
nitude of 20.53. The panels in Figure [1] show from left 
to right cut-outs from a pre-event image (2001 Novem- 
ber 20), an image 4 days after discovery (~6 days before 
the peak; 2003 December 17), and their difference image, 
respectively. The event position in the difference images 
is 0.064 ± 0.012 arcsec to the south of a source (the pre- 
sumed host galaxy) in the pre-event template image (see 
Figured]). At a redshift of z = 0.289 of SN 2003ma 
and its host (see Section [2~3|l . this corresponds to a pro- 
jected distance of d m 460 ± 85 pc. The host galaxy 
has apparent magnitudes of Bsm = 20.93 ± 0.06 mag, 
VRsm = 20.15 ±0.02 mag, and I SM = 19.76±0.03 mag. 

SN 2003ma was also detected by the OGLE-III sur- 
vey (2001 - 2009), which is a photometric survey using 
a dedicated 1.3-m Warsaw telescope located at the Las 
Campanas Observatory, Chile, operated by the Carnegie 
Institution of Washington. The camera uses eight SITc 
2048 x 4096 CCD detectors with 15 fim pixels resulting in 
0.26 arcsec pixel -1 scale and 35 x 35 arcmins total field 
of view. For the details of the instrumentation setup 



< 



0.15 - 



0.10 - 



0.05 



0.00 



-0.05 



• Difference Images 
■ Pre-event Images 




-o.io 



-0.05 



o.oo 

ARA(") 



0.05 



0.10 



Fig. 2. — Centroids of the sources at /near the position of 
SN 2003ma from multiple images. The red squares indicate the 
centroids of detections in the pre-event (MJD < 52965) Super- 
MACHO images, while the blue circles indicate the high S/N cen- 
troids of SN 2003ma in the SuperMACHO difference images during 
52985 < MJD < 53020. SN 2003ma is offset from the template 
source by 0.064 ± 0.012 arcsec, indicating that it is at a projected 
distance of d fs 460 ± 85 pc from the host. 



we refer the reader to lUdalskil (|2003t ). Approximately 
500 photometric points per star towards the LMC were 
accumulated over eight seasons, between July 2001 and 
May 2009. Most of the observations were carried through 
the Cousins /-band filter with exposure times 180 sec- 
onds. Full description of the reduction techniques, differ- 
ence image analysis, photometric calib ration and astro- 
metric transformations can be found in lWozniakl (|2000D : 
lUdalski et all (12001) . 

Pre-event /-band data were also gathered during the 
OGLE-II survey (1996 - 2000), which used the same 
dedicated telescope as OG LE-III, but was equ ipped with 
a single 2048 x 2048 CCD (|Udalski et alJfl997h . Since the 
QE of the two CCD's are slightly different, we apply a 
small 2% correction to the OGLE-II data so that their 
average flux matches the pre-event flux in the OGLE-III. 
Since the OGLE-II data is several years before the event, 
we used it to test for pre-event variability, for example 
caused by an active galactic nucleus (AGN). 

The upper and lower panels of Figure [3] shows the 
difference image light curves from the SuperMACHO 
and OGLE surveys. Images with poor seeing and/or 
significant cloud cover are more likely to have large 
(and non-Poissonian) systematic uncertainties. There- 
fore, only images with seeing better than 1.5" and with 
atmospheric extinction less than 1 mag are used for the 
analysis. We make an exception for images from July 
and August 2008, which typically have poor seeing, but 
are crucial for understanding the late-time behavior of 
SN 2003ma. For these months, images with seeing better 
than 1.9" were retained. Since there are several images 
in all filters from 2008, we were able to average the flux 
from these images. 
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Fig. 3. — Difference image light curves of SN 2003ma. Top: 
Black triangles, blue circles, and orange diamonds show the 
SuperMACHO difference flux light curves in Bsm, VHsMi an( i IsMi 
respectively, from 2001 November 11 to 2008 August 25. Middle: 
Small black circles indicate the Tqgle difference flux light curves. 
Bottom: Because of the low S/N of a single detection, we average 
the data (orange squares). The horizontal dashed lines indicate 
zero flux. The vertical dotted line indicates the transition between 
OGLE-II and OGLE-III measurements. 



Because difference imaging reveals more instrumental 
and reduction artifacts than there are real variable ob- 
jects in an image, st andard profile-fitting s oftware pack- 
age, like DoPHOT (jSchechter et all 119931) . have prob- 
lems determining the optimal PSF in the difference im- 
age. To remedy this, we implement a customized ver- 
sion of DoPHO T that uses the applicable PSF fr om the 
original image (jRest et al.ll2005tlGarg et al.H2007t) . After 
an object is detected in any difference image we "force" 
photometry in all difference images at the event posi- 
tion regardless of the flux level. We thus measure the 
flux in pre-event images and negative fluxes are possi- 
ble. These measurements are critical for constraining 
the light curves. For example the pre-event light curve 
(MJD < 52970) is consistent with no variability, indi- 
cating that the host galaxy has not demonstrated a his- 
tory of variability that would indicate the presence of an 
AGN. 

In addition to the SuperMACHO VR images, several 
B and / images were obtained on 2001 November 20, 
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Fig. 4. — Optical light curves of SN 2003ma. The black trian- 
gles, blue circles, orange diamonds, and orange squares show the 
apparent magnitudes in the natural system of -Bsm> V-Rsm> ^SM> 
and /qglEi respectively. Because of the low S/N of a single detec- 
tion, we average /ogle (orange squares), with the range of epochs 
used to calculate this average indicated by the error bars along 
the time axis. The inset shows the light curve during peak from 
52928 < MJD < 53130, and the 3a brightness upper limits. 



2004 September 20, 2004 October 21, 2008 July 27, and 
2008 August 26. The images from 2001 November 20, ob- 
tained before the event, are used as the template images 
for these filters. 

The middle panel of Figure |3] shows the OGLE I- 
band difference image photometry (small black circles). 
Mainly because of the difference in telescope aperture 
(1.3 m for OGLE and 4 m for SuperMACHO), the signal- 
to-noise ratio (S/N) of the OGLE photometry is sig- 
nificantly lower than that of the SuperMACHO data. 
To obtain measurements with sufficient S/Ns, we have 
binned the OGLE data (lower panel of Figure [3]). Since 
the OGLE observations are taken in blocks, we divide 
each block in two equal parts, each containing the same 
amount of measurements. The only exception is dur- 
ing peak, where the S/N of the measurements is signif- 
icantly higher, and in addition the light curve evolves 
faster. Thus around peak (52975 <= MJD <= 53015), 
we average the data in 5-day bins, and in 25-day bins 
in the range 53015 < MJD <= 53150. In the following 
analysis, we only use the averaged OGLE data. 

We convert the CTIO 4 m difference flux detections 
into apparent magnitudes (see Figure |4] and Table [T]) in 
the natural CTIO 4 m magnitude system, and we correct 
for extinction as described in Section 14.11 The OGLE 
measurements have been converted into standard I-band 
magnitu des using the s t andard OGLE phot ometric cali- 
bration (|WozniaJdl20txi lUdalski et al.ll2008f ). 

2.2. Infrared Photometry 

We use archived images and catalogs from the 
Spitzer Space Teles cope SAGE survey of the LMC 
(jMeixner et al.l 120061 ) to measure the mid-infrared pho- 
tometry of SN 2003ma. The SAGE survey covered a 
7° x 7° region centered on the LMC during two epochs, 

2005 July 15 - 26 (epoch 1) and 2005 October 26 : 
November 2 (epoch 2) for its IRAC observations, and 
2005 July 27 - August 3 (epoch 1) and 2005 November 



The Extremely Energetic SN 2003ma 



TABLE 1 

Apparent magnitudes of SN 2003MA 



MJD 



VRsi 



•foGL 



52959 
52964 
52966 
52970 
52972 
52974 
52976 
52979 
52982 
52985 
52986 
52988 
52988 
52990 
52992 
52992 
52994 
52996 
52997 
52998 
53000 
53001 
53004 
53011 
53019 
53029 
53052 
53077 
53106 
53257 
53268 
53268 
53285 
53289 
53295 
53297 
53299 
53299 
53299 
53315 
53325 
53329 
53331 
53348 
53350 
53352 
53354 
53356 
53404 
53414 
53639 
53641 
53649 
53652 
53654 
53663 
53684 
53704 
53706 
53712 
53714 
53733 
53738 
53777 
54006 
54127 
54402 
54502 
54674 
54674 
54704 
54704 



29126 
27804 
28245 
30195 
29799 
25935 
26257 
64149 
71623 
24459 
27416 
68744 
28817 
25875 
29069 
73489 
31115 
28857 
44242 
25951 
31187 
72660 
32118 
28573 
31815 
29812 
97253 
45696 
03597 
38920 
26639 
27047 
33774 
21636 
33515 
22518 
36194 
36544 
36742 
28817 
24440 
29664 
34066 
18819 
32545 
23171 
23758 
25495 
63460 
22055 
30253 
33128 
38498 
32054 
22569 
32678 
22316 
30592 
22709 
35217 
35882 
16559 
16874 
61556 
82171 
67356 
73904 
13612 
41064 
42468 
35730 
37563 



23.15 



23.34 



(0.17) 



(0.40) 



< 24.19 

< 24.01 

< 24.16 



24.22 
23.46 
22.50 



20.95 
20.81 

20.24 
19.83 
19.82 

19.81 
19.87 

19.95 
20.03 

20.12 

20.23 



0.43) 
0.19) 
0.09) 



0.03) 
0.03) 

0.02) 
0.02) 
0.03) 

0.02) 
0.02) 

0.02) 
0.02) 

0.02) 

0.02) 



22.02 (0.05) 



22.22 
22.08 
22.18 
22.15 



22.08 
22.12 
22.17 
22.11 
22.20 
22.17 
22.21 
22.18 
22.16 

22.11 
22.77 
22.60 
22.68 
22.81 
22.68 

22.67 
22.65 
22.56 
22.69 
22.69 
22.62 
22.71 



0.06) 
0.06) 
0.06) 
0.06) 



0.06) 
0.06) 
0.06) 
0.06) 
0.06) 
0.07) 
0.06) 
0.06) 
0.07) 

0.06) 
0.10) 
0.12) 
0.11) 
0.11) 
0.10) 

0.09) 
0.09) 
0.10) 
0.12) 
0.29) 
0.08) 
0.09) 



23.82 (0.41) 
24.19 (0.20) 



21.15 



21.13 



(0.06) 



< 22.30 



< 21.05 

< 19.80 



20.21 

19.65 

19.77 

19.83 
19.92 



(0.18) 

(0.10) 

(0.09) 

(0.08) 
(0.10) 



19.90 (0.04) 

19.91 (0.05) 
20.02 (0.05) 
20.53 (0.07) 



21.05 



(0.10) 



21.30 



21.24 



(0.07) 



(0.12) 



(0.08) 



21.39 (0.10) 
21.87 (0.20) 

21.46 (0.14) 
22.26 (0.22) 

22.47 (0.23) 



22.61 (0.45) 
22.59 (0.15) 
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2-9 (epoch 2) for its MIPS observations. We used the 
epoch 1 and 2 images provided by the SAGE team to 
produce small (300 pixel wide) cutouts around the po- 
sition of SN 2003ma, for all four IRAC bands and the 
24 urn MIPS band. The IRAC and MIPS cutouts have 
spatial scales of 1.2 and 2.5 arcsec pixel -1 , respectively, 
and are all on a common coordinate reference frame. 

Although an object was clearly visible at the position 
of SN 2003ma on several of the Spitzer images, we found 
that the 5^4 GE catalog did not contain photometry for it 
in all bands, likely due to the presence of several bright 
stars near the SN 2003ma position. Moreover, because 
the crowding in the Spitzer images forced us to conduct 
artificial star tests in order to understand the true photo- 
metric errors, we performed our own PSF-fitting photom- 
etry on the im age cutouts using DAOPHOT / ALLSTAR 
(|Stetsonlll987t) . We used a threshold of ltr to find stars 
and a sequence of apertures for initial photometry. To de- 
rive a PSF for each image, we first selected 30 - 130 PSF 
stars for the IRAC images, and 18 for the MIPS 24 fj,m 
image. We then used an iterative procedure of fitting a 
PSF function to the PSF stars, subtracting neighboring 
sources, and then refitting the PSF function to the PSF 
stars, while increasing the order of the spatial variability 
of the function up to a maximum of two. After several 
iterations, we used the derived PSF to derive photome- 
try for all detected sources. After subtracting all sources 
except for the PSF stars, we calculated aperture photom- 
etry of the PSF stars out to a radius of 11 pixels (13'.'2) 
for the IRAC images, and 29 pixels (7275) for the MIPS 
24 /im images. We used this aperture photometry to ap- 
ply a correction to the PSF-fitted photometry. Finally, 
we matched the aperture-corrected photometry to that of 
sources in the SAGE catalog, and applied the zero-point 
differences to provide the photometric calibration. 

For three of the images (epochs 1 and 2 for the 3.6 /jm 
band, and epoch 1 for the 4.5 /j,m band), we found that 
the crowding in the images and the intrinsic faintness 
of the object left SN 2003ma still unmeasured, with no 
matching sources within 3 pixels of its position. For these 
three images, we thus use the 8.0 /zm epoch 1 image 
positions as the master source list, and forced ALLSTAR 
to fit photometry at the positions, allowing only the sky 
and object magnitude to be free parameters. 

To ensure that our derived photometric uncertainties 
were reasonable in the Spitzer images which all suffer 
from crowding, we performed artificial star tests. Arti- 
ficial stars were added, a few at a time, to copies of the 
original images, measured using the photometric proce- 
dure described above, and then analyzed to evaluate the 
photometric errors, offsets, and completeness as a func- 
tion of input magnitude. In the course of these tests, 
we added a total of 10,000 artificial stars, with magni- 
tudes distributed evenly over the range 14 < M 3 . 6 Mm < 
18 mag, and with colors with respect to the other bands 
matched to those observed for SN 2003ma. The stars 
were added 50 at a time (—5% of the real detected sources 
in each image) at random positions. Following the pho- 
tometric measurement step, we recovered the detected 
artificial stars by first removing from the output lists all 
stars matching real sources detected in the original im- 
ages, and then comparing the remaining sources with the 
list of input artificial stars. We then computed the av- 
erage offset and dispersion in the input minus recovered 



magnitudes and the fraction of artificial stars recovered 
as a function of input magnitude. 

Tablc[2]lists the Spitzer photometry of SN 2003ma, and 
the quantities derived from the artificial star tests. The 
measured flux in mJy is shown in the 3rd column. The 
flux uncertainty in the 4th column was determined by 
DAOPHOT and is typically between 5% - 15%. This is 
the appropriate relative uncertainty between the fluxes 
measured in the same band but different epochs, since 
any bias due to crowding should have nearly the same 
effect in both epochs assuming that the neighboring ob- 
jects are not variable. For different wavelengths, the 
crowding is different and thus has a different impact 
on the photometry. Therefore our estimate for the to- 
tal photometric uncertainty (6th column in Table [2]) of 
SN 2003ma at each epoch is based on the mean photo- 
metric uncertainty arrived at in the appropriate magni- 
tude range of the artificial stars lists. For all comparisons 
of fluxes from different bands, we use these total errors, 
which are -30%, -30%, -30 - 50%, -10 - 15%, and 
-20% for the 3.6, 4.5, 5.8, 8.0, and 24 urn images, re- 
spectively. 

The average absolute offsets in the input minus re- 
covered magnitudes of the artificial stars are typically 
smaller than 0.1 mags, such that we do not expect the 
Spitzer detections of SN2003ma to be spurious. The only 
exception is the 5.8 micron photometry, for which the ab- 
solute offsets are 0.5 mags. The average completeness for 
the 5.8 /im images at the magnitude of SN 2003ma was 
also <50%, which suggests that the 5.8 fxm photometry 
is unreliable compared to that in the other bands. As 
a consistency check, we compared our unforced photom- 
etry to the recently complete, final catalog and ar chive 
lists from the SAGE project (jMeixner et al.l 120061 ) and 
found them to be in agreement. 

2.3. Spectroscopy 

We followed SN 2003ma spectroscopically for nearly 
six year s using the 8 m Gemini-South telescope (with 
CMOS; IHook et alJl200l and the 6.5 m Magellan Clay 
telescope (with LDS S2; lAllington-Smith et al.l I1994L 
LDSSSE3, and MagE: IMarshall et al.1 l200l . The spec- 
traFI are presented in Figure [5] and a journal of observa- 
tions are presented in Table |3] 

Standard CCD processing and spectrum extraction 
were accomplished with IRAF. The data we re extracted 
using the optimal algorithm of iHornel (|1986f ). Low-order 
polynomial fits to calibration-lamp spectra were used to 
establish the wavelength scale, and small adjustments de- 
rived from night-sky lines in the object frames were ap- 
plied. For the MagE spectrum, the sky was sub tracted 
from t he images using the method described by iKelsonl 
(|200l . We employed our own IDL routines to flux cal- 
ibrate the data and remove telluric lines using the well- 
exposed continua of the spectrophotometric standards 
dWade fc Hornd[l988l iFblev et al.ll2003l I2009T) . 

The 2004 spectrum was obtained far from the parallac- 
tic angle at high airmass. As a result, its continuum and 

30 http://www.lco.cl/telescopcs- 
information / magcllan / instruments- 1 /ldss-3- 1 / . 

31 The Gemini South GMOS spectra were obtained as part of 
the GS-2003B-Q-12 science program and the header target name 
in the images and spectra was candl0194.sm76_4. 
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TABLE 2 

IRAC and MIPS Photometry and Artificial Star Test Results 



Band 


t a 


Flux b 


Flux relerr c 


AFlux d 


Flux err e 


— 

Compl. 


(|im) 




(mjy) 


fin Iv) 


Cm Ivl 


im Ivl 


(%) 


3.6 


1 


0.063 g 


0.006 


0.006 


0.018 


46.4 


1.5 


1 


0.048s 


0.003 


0.003 


0.014 


53.3 


5.8 


1 


0.059 


0.004 


0.017 


0.017 


38.9 


8.0 


1 


0.108 


0.004 


0.014 


0.018 


77.3 


24.0 


1 


1.164 


0.080 


0.099 


0.227 


94.4 


3.6 


2 


0.055 s 


0.003 


0.005 


0.019 


40.9 


1.5 


2 


0.065 


0.003 


0.006 


0.014 


65.6 


5.8 


2 


0.028 


0.004 


0.015 


0.013 


13.6 


8.0 


2 


0.184 


0.006 


0.014 


0.022 


89.3 


24.0 


2 


1.141 


0.108 


0.098 


0.228 


94.0 



a Epoch ID. Epoch 1 spans from 2005 July 15 - August 3, and epoch 2 from 2005 
October 26 - November 9. 

k Flux at the position of SN 2003ma. This flux contains both the flux of SN 2003ma 
as well as any potential flux from the host galaxy. 

c Flux error as determined by DAOPHOT. This is the appropriate relative error 
between the fluxes measured in the same band at different epochs 
^ Average detected flux minus input flux as derived from artificial star tests, inter- 
polated at the measured flux level of SN 2003ma in each image. 

e Average dispersion in the detected flux minus the input flux from the artificial 
star tests, interpolated at the measured flux level of SN 2003ma in each image. 
^ Average completeness in %. 
s Forced photometry. 
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Fig. 5. — Optical spectra of SN 2003ma. The spectra are de- 
noted by the year in which they were obtained. The spectra in 
years 2003, 2004, 2005, and 2008 correspond to rest-frame phases 
of 1, 279, 905, and 1423 days relative to B maximum, respectively. 
There is significant galaxy contamination in all spectra, likely dom- 
inating the spectrum at later phases. Since the 2004 spectrum was 
not obtained at the parallactic angle, its spectral shape is likely 
incorrect. 



TABLE 3 

Log of Spectral Observations 



Phase a 


UT Date 


Telescope / 


Exposure 






Instrument 


to 


1 


2003 December 22.3 


Gemini South/GMOS 


3 x 610 


279 


2004 December 14.1 


Magellan/LDSS2 


1500 


905 


2005 October 29.3 


Magellan/LDSS3 


900 


1423 


2008 December 28.3 


Magellan/MagE 


2 x 1800 


1649 


2009 October 17.3 


Magellan/MagE 


1800, 1200 



a Rest-frame days since B maximum, 2003 December 20.8 (JD 2,452,994.3). 

any derived line ratios are likely incorrect (|Filippenkol 



Fig. 6.— Optical spectrum of SN 2003ma from 2008 (i = 
1423 days). Narrow emission lines are marked. The continuum 
is likely dominated by the host galaxy, while the broad Ha and 
potential He II A4686 lines are from SN 2003ma. 



[19821) . Since the offset between the SN and host is very 
small, we were not able to perform local galaxy subtrac- 
tion during extraction. For the 2003 and 2008 spectra, 
we have calibrated our absolute spectrophotometry to 
match concurrent photometry of the SN and host com- 
bined. This is not possible for the 2004 and 2005 spectra 
because of incorrect relative spectrophotometry and a 
small wavelength range that does not fully cover any of 
our broad-band filters, respectively. 

Several narrow galactic emission lines are present in all 
spectra and are identified in Figure [5] at a redshift of z — 
0.289. Throughout this paper we assume a cosmology 
with H = 70 kms" 1 Mpc" 1 , D, m = 0.3, tt A = 0.7, 
which yields a luminosity distance of c?l = I486 Mpc for 
the SN redshift. 
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3. COMPARISONS TO OTHER CORE-COLLAPSE 
SUPERNOVAE 

3.1. Spectroscopic Comparisons 

The five spectra of SN 2003ma are presented in Sec- 
tion 12.31 The 2008 and 2009 spectra are dominated 
by host-galaxy light (see Figure E]); however, there is a 
strong, broad component to the Ha line. We are confi- 
dent that this broad component is from the underlying 
SN and not from another source such as an AGN (for 
instance, the broad component is not visible in our 2003 
spectrum). 

The 2003 spectrum is much bluer than the late-time 
spectra, indicating that it contains more SN light than 
the late-time spectra, which is not surprising consider- 
ing that the SN was ^3 mag brighter in 2003 compared 
to 2008. In order to isolate the SN continuum in the 
2003 spectrum, we have subtracted a galaxy template 
spectrum from the 2003 spectrum. We chose the galaxy 
template to match the general continuum shape of the 
late-time spectra. The pre-event photometry can be used 
to properly scale the galaxy template. Similarly the pho- 
tometry at the time of each SN spectrum can be used 
to scale each spectrum. Unfortunately, this procedure 
produces a galaxy spectrum that is ~ 10% too bright 
compared to the late-time spectrum (assuming that the 
late-time spectrum is dominated by the host galaxy). 
However, considering the fact that the galaxy template 
spectrum does not include emission lines and there may 
be crowding in the template image, we do not consider 
this to be a significant discrepancy. The galaxy template 
spectrum was rescaled to match the continuum of the 
2008 spectrum and subtracted from both the 2003 and 
2008 spectra. 

We present the dereddened (our spectra of SN 2003ma 
have only been corrected for foreground extinc tion based 
on Galactic dust maps ISchlegel et aD Il99l see Sec- 
tion l4.1j) galaxy-subtracted 2003 and 2008 spectra in Fig- 
ure [TJ We compare these spectra to those of SN 2005ip 
(|Smith et al.l l2009bh . after dereddening the spectra of 
SN 2005ip by E(B - V) = 0.047 mag to account for 
Milky Way extinction. In order to match the contin- 
uum slope of SN 2003ma, we add an additiona l E(B — 
V) = 0.9 mag of host extinction. iSmith et all (|2009bl ) 
found that SN 2005ip may have had host extinction of 
E(B — V) = 0.3 mag based on measurements of Na D 
absorption, however since the absorption features were 
only clearly detected in the spectrum taken 1 day after 
discovery, they were unsure of the true magnitude of ex- 
tinction. It seems unlikely that SN 2005ip suffered from 
E(B — V) = 0.9 mag of reddening without strong Na D 
absorption. Instead, the difference in continuum shapes 
may indicate that SN 2003ma was hotter at this epoch 
than SN 2005ip (see also Section 143]) . 

At maximum, both SN 2003ma and SN 2005ip have 
relatively featureless continua. With the applied extinc- 
tion corrections, there are no significant deviations be- 
tween the two spectra. The spectrum of SN 2003ma 
is consistent with that of a blackbody. There is a low 
amplitude, broad feature in the SN 2003ma spectrum 
at ^6150 A(in the rest frame) that may correspond to 
Ha absorption. If this identification is correct, the fea- 
ture would have an absorption minimum blucshiftcd by 
^20,000 kms -1 . In Figure [3 we also compare the 2008 
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Fig. 7. — Optical spectra of SN 2003ma compared with those of 
SN 2005ip. Phases are given in rest-frame days since maximum 
light. Top: Galaxy-template subtracted maximum-light spectrum 
of SN 2003ma (black) and maximum-light spectrum of SN 2005ip 
(blue). An additional E(B-V)=0.9 mag of internal extinction was 
applied to SN 2005ip to match the continuum shape of SN 2003ma. 
Bottom: The 1423 day (4 year; black) and 905 day (2.5 year; blue) 
spectra of SN 2003ma and SN 2005ip, respectively. A template 
galaxy spectrum has been added to the SN 2005ip spectrum to 
match the continuum shape of the SN 2003ma spectrum. 



spectrum of S N 2003ma to the t = 905 day spectrum 
of SN 2005ip ([Smith et al.ll2009bO . The continua of the 
spectra match, but this may be the result of an incorrect 
treatment of the galaxy contamination. The only feature 
of significance is the Ha line, which is much broader for 
SN 2003ma compared to SN 2005ip. The broad Ha line 
is also seen in the 2009 spectrum, but at lower equivalent 
width. 

We detect narrow He I A5876 emission in the spectra 
from 2004, 2008, and 2009, but do not detect it in the 
2003 spectrum (the wavelength range of the 2005 spec- 
trum does not cover this feature) . The 2003 spectrum has 
a relatively low S/N, and the He I may be too weak to be 
detected. Alternatively, in Section 14*^21 we argue that the 
blackbody temperature in 2003 was > 15,000 K, which 
should produce sufficient high-energy photons to singly 
ionize He. In Figure |6l we tentatively identify relatively 
broad He II A4686 in the 2008 spectrum. If correct, this 
would require a significant amount of ionizing far-UV or 
X-ray photons from the circumstcllar interaction to ion- 
ize He in the post-shock material. However, as shown 
in Figure [7j SN 2005ip has several narrow lines in this 
region of the spectrum. It is possible that the 2008 spec- 
trum of SN 2003ma does not have a high enough S/N to 
see these individual lines. 

In SNe Iln, the Ha line can have up to four distinct 
components : broad , intermediate, narrow, and galactic. 
ISmith et all (|2009bD was able to identify the first three 
regions for SN 2005ip, while the fourth component was 
not seen. Physically, the broad, intermediate, and nar- 
row components are formed in the fast-moving SN ejecta, 
the post-shock circumstcllar material, and the pre-shock 
circumstellar material, respectively. The widths of the 
broad and narrow components indicate the velocity of 
the SN ejecta and the circumstellar wind/ejection, re- 
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spectively, while the intermediate component is an indi- 
cation of the amount of energy transferred from the shock 
to the circumstellar material. The galactic component, 
which is the line-of-sight emission from the host galaxy 
and physically unrelated to the SN, may have a similar 
velocity to that of the narrow component . 

In Figure [5] we show the narrow [O I] A6300, Ha, 
[N II] AA6548, 6583, and [S II] AA6716, 6731 emission 
lines fitted with a Voigt profile with its width fixed to 
that of the [N II] A6583 line in each spectrum. The Voigt 
profile produces better fits to the wings of the narrow 
lines than a single Gaussian. The narrow lines are unre- 
solved in the 2003 and 2004 spectra, but the high reso- 
lution of the 2008 and 2009 spectra resolves the narrow 
lines to have a FWHM = 110 kms" 1 . In 2004, 2008, 
and 2009 we fit the broad base beneath the Ha+[N II] 
narrow-line complex with a Gaussian. The best fit was 
for a Gaussian with FWHM = 5800 ± 800 kms" 1 and 
a velocity offset of +350 kms -1 relative to the narrow 
Ha line. This velocity offset could be due to asymme- 
tries in the line profile. The broad component is not 
detected in the 2003 spectrum, and we place an upper 
limit on the presence of this Gaussian component to be 
< 0.2 times the flux of the narrow Ha line. The ratio of 
the broad component to the narrow component of Ha is 
0.7 ±0.1, 0.4 ±0.2, and 0.5 ±0.2 in 2004, 2008, and 2009, 
respectively The change in this ratio is not statistically 
significant, and systematic effects such as different posi- 
tion angles of the slit and varying seeing can also produce 
differences. 

We measure the standard diagnostic [O III]/H/3 
and ]N II] /Ha, and JS II] /Ha narrow- line ra - 
tios (jBaldwin et all 119811: IVeilleux fc Osterbrocklll987t ). 
which are insensitive to extinction and atmospheric re- 
fraction, for each spectrum from the Voigt profile fits. 
The narrow-line ratios, listed in Table IH are all con- 
sistent with a star-forming galaxy. This, in combina- 
tion with the lack of variability in the densely sampled 
SuperMACHO and OGLE difference imaging data over a 
baseline of 7 years before the SN, place strong constraints 
on the presence of an AGN in the host galaxy. 

If the broad Ha component is powered by CSM inter- 
action, then we expect there to be some contribution to 
the narrow Ha and H/3 emission lines. However, with 
no observed velocity offset between the expected narrow 
component from the SN and the host galaxy emission 
lines, it is difficult to tell if the SN contributes to the 
narrow emission lines. Furthermore, the narrow-line ra- 
tios are measured to be constant within the errors from 
2003 to 2008, and thus the lines are most likely domi- 
nated by the H II regions in the host galaxy. The lumi- 
nosity of the Ha line in 2003, when the broad compo- 
nent of Ha associated with the SN is not yet observed, 
is L(Ha) = 6.1 x 10 42 ergs s -1 , which corresponds to a 
substantial star- formation rat e of 4.8M,?) year -1 charac- 
teristic of a starburst galaxy (jKennicuttl 1 1 9981 ) . 

At day 905, SN 2005ip had been on its plateau for 
over 2 years, and the broad component of its Ha line had 
faded by a factor of ~6 while its narrow -(-intermediate Ha 
line, which is assumed to be powered by interaction with 
the CSM, had increased by a facto r of 3, resulting in a 
nearly constant total Ha luminosity (jSmith et al.ll2009bl ). 
At this epoch, the FWHM of the intermediate component 



and broad component of Ha wer e 740 ± 250 km s 1 and 
8800 ±660 kms" 1 , respectively (jSmith et al.ll2009bft . In 
comparison, at day 279 and 1423, SN 2003ma was on its 
second plateau and also most likely dominated by cir- 
cumstellar interaction, however the width of the broad 
component of its Ha line is at least 5 times larger than 
the intermediate component seen in SN 2005ip. Sim- 
ilarly, it is significantly broader than the intermediate 
component of Ha observed in SNe Iln 1988Z, 2006tf, and 
2007rt which had a characteristic width of 2000 kms -1 
(lAretxaea et al.lll999l : iSmith et aLll2008at iTrundle et all 
l200l " and SN 2006gy w hich had an interme diate compo- 
nent with 4000 kms" 1 (|Smith et al.ll2007b| ). 

The ~6000 kms -1 component seen in the 2004 and 
2008 spectra of SN 2003ma could be associated with SN 
ejecta that has decelerated over time. The width of the 
broad component of Ha was observed to decrease from 
~15,000 kms " 1 on day 34 to 5000 k ms" 1 after day 500 
in SN 1988Z (lAretxaga et al.l 11999ft and decrease from 
~10,000 kms -1 o n day 102 to < 6000 kms" 1 after day 
273 in SN 2007rt (|Trundle et al.l 120091 ) . It may be that 
the day 1 spectrum of SN 2003ma was too early, and 
the day 279 spectrum was too late to catch the broad 
component when it had a velocity width of tens of thou- 
sands of kms" 1 . However, given the long-lived plateau 
in the light curve, circumstellar interaction must be the 
dominant source of energy after day 90. Thus, the broad 
component in SN 2003ma is most likely an exceptionally 
fast "intermediate component" that is the result of the 
interaction with the CSM of a highly energetic explo- 
sion, even more energetic than SN 2006gy (see Section 
14. 5[) . which imparts a much faster blast-wave speed to 
the postshock shell. However, electron scattering can 
also be a significant source of broadening of the Ha line 
in interacting SNe, and thus the speed of the material 
may be slower than inferred from the width of the line 
(jDessart et al.ll2009f ). 

3.2. Photometric Comparisons 

The I-band light curve of SN 2003ma has two nearly 
flat plateaus, and strong emission in the V and I bands 
for ~ 1000 days, reminiscent of the long-lived emission 
seen in interacting SNe Iln. Although the first plateau 
in the rest-frame R-band is similar to the initial plateau 
seen in SNe IIP, its comparable plateau in the rest- 
frame B band is unlike IIP SNe which begin their de- 
cline immediately after the peak in the B band. Fig- 
ure [TT] shows a c omparison of SN 2003ma to canonical 
SN I IP 1999em (|Leonard et al.l 120021: lElmhamdi et al.l 
I2003D. along with ex treme interac ting SNe Iln 1988 Z 
(|Aretxaga etaIl[T999h and 2005ip (ISmith et al.ll2009bl) . 
and extremely luminou s SN Iln 2006g y (IQfek et al 



2007[ ISmith et al.ll27)07bh. SN IIL 200 5ap ( Quimbv et al 
20071) SN IIL 2008es dMiller et al.l [20091: iGezari et al 



20093), and normal SN la 1988aq (|Riess et al.ll2005D . All 



light curves are corrected for Galactic extinction, and 
in the case of SN 1999em and SN 2006gy for int ernal 
extinction (jLeonard et al.ll2002t ISmith et al.H2007bD . K- 
corrections are applied to SNe at distances farther than 
100 Mpc. 

In SNe IIP, a plateau is observed in the V and R bands 
as a result of H recombination in the SN ejecta. The 
plateau is not seen at shorter wavelengths (U and B 
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TABLE 4 
Diagnostic Narrow-Line Ratios 



Year 


[OIIIJA5007/H/3 


[NII]A6583/Ha 


[SII]AA6716, 6731/Ha 


2003 


1.98 ± 0.07 


0.18 ± 0.01 


0.21 ± 0.02 


2004 


2.19 ± 0.03 


0.23 ± 0.04 


0.21 ± 0.03 


2005 


2.10 ± 0.07 






2008 


2.19 ± 0.05 


0.22 ± 0.02 


0.18 ± 0.01 


2009 


1.71 ± 0.12 


0.19 ± 0.02 


0.17 ± 0.01 
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Fig. 8.— Narrow [O I] A6300, Ho, [N II] AA6548, 6583, and [S II] 
AA6716, 6731 emission lines fitted with Voigt profiles with a width 
equal to the [N II] A6583 line in each spectrum, respectively, are 
shown with a blue line. In 2003 the continuum if fitted with a 
quadratic function, and in 2004, 2008, and 2009 the continuum is 
fitted with a constant value. An additional broad Gaussian com- 
ponent that is needed to fit the Ha line in 2004, 2008, and 2009 is 
shown in red. In 2004 the [N II] A6548 line is not detected above 
the noise, and is weaker than expected for a line ratio of 3:1 for 
[N II] A6583/A6548. 



bands) due to sensitivity to the declining temperature 
and line blanketing from metals, nor at longer wave- 
lengths (7 band) since the bandpass is o n the Rayleigh- 
Jeans tail of the emission the whole time ([Leonard et al.1 
l200l . 

The SuperMACHO data constrains the pre-maximum 
light curve of SN 2003ma exceedingly well (see Fig- 
ure [T2|) . The first significant detection is 4.5 mag be- 
low peak 20 days before maximum, and we estimate 
the day of explosion to be between the last upper limit 
and the first detection at t = MJD 52,971 ± 1. The 
slow decline of the I band (rest-frame R band) in the 
first plateau (53010 < MJD < 53110 or t ~ t = 
40 - 140 days) of 0.0017 ± 0.0009 mag day" 1 and in the 
second plateau (MJD > 53110 or t-t = 140-1735 days) 
of 0.0011±0.0001 mag day" 1 is much flatter than the 481 
day decline seen in S N Iln 2007rt with 0.006 mag day -1 
(|Trundle et all 120091 ), and is more like the extreme 
plateaus observed in the light curves of SNc Iln 1988Z 



and 2005ip (see Figure fT2l . The slow decline in these 
SNc Iln light curves are explained by an excess of emis- 
sion above the photospheric emission produced by the 
interaction of the SN ejecta with circumstellar material. 
The spectra of these SNe demonstrated broad Balmer 
emission lines, with a narrow component and no P-Cygni 
absorption, suggesting that the emission is confined to an 
outer shell (jTuratto et al.lll993D . 

The color evolution of SN 1988Z and SN 2003ma, 
shown in Figure 1131 arc remarkably similar. The ex- 
tremely flat plateaus of this type of SN are attributed 
to an increase in narrow and intermediate- width Ha flux 
as well as a forest of othe r narrow lines powe red by the 
circumstellar interaction (jSmith et al.ll2008al) . The de- 
creasing Ha line velocity wid th over time in SN 19 88Z 
was successfully modeled by lAretxaga et al.l (|1999|) as 
the consequence of a SN remnant shock expanding into 
a dense medium. The fact that SN 2003ma is bluer than 
SN 1988Z in particular during the early phases can be 
attributed to the higher blackbody temperature. 

4. PHYSICAL PARAMETERS OF SN 2003ma 

4.1. Extinction 

We use the typ ical galactic extinction law of 
ICardelli et al.1 (|1989f ) parameterized by a value of Ry = 
A V /E(B - V) = 3.1. We obtaii£3 a combined Galactic 
and LMC reddening of E( B — V) = 0.348 mag at the 
position of SN 2003ma from lSchlegel et al.l (|1998[) . As a 
consistency check, we examine the spectra of SN 2003ma 
for tracers of interstellar dust. The Na D doublet is not 
detected at or near zero velocity in any of our spectra. 
The 2008 spectrum of SN 2003ma, which has a resolu- 
tion of ~20 kms -1 , is ideal for examining the Na D dou- 
blet. We present the wavelength regions from that spec- 
trum which correspond to the rest-frame Na D doublet 
for the Milky Way, LMC, and host galaxy of SN 2003ma 
in Figure [9] From this spectrum, we can rule out Na D 
absorption with an equivalent width of > 0.4 A from 
the Milky Way and the LMC. Using the relationship of 
iBarbon et all (| 19901 ). we expect E(B - V) < 0.1 mag 
from the Milky Way and the LMC, respectively. How- 
ever ,|BlondEeEaD|200i) found E(B-V) > 0.3 mag for 
objects with Na D equivalent widths below about 0.3 A. 
We therefore believ e that our mea s urem ents of Na D are 
consistent with the lSchlegel et al.l (|1998t ) extinction esti- 
mates. 

In Figure [9l we see that there is a possible detection of 
Na D absorption from the host galaxy at a redshifted ve- 
locity of v = 200 kms -1 relative to the host galaxy with 
a width of < 20 kms -1 , the resolution of our spectrum. 
However, we do not see a narrow emission line compo- 
nent from the SN that is offset from the emission lines of 

32 http:/ /irsa. ipac.caltech.edu/applications/DUST/ 
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Fig. 9. — Optical spectrum of the 1423 day spectrum of 
SN 2003ma near the rest wavelength of the Na D doublet for the 
Milky Way (top), LMC (middle), and host galaxy (bottom), re- 
spectively. The vertical dashed lines indicated the wavelengths of 
the Na D lines. We detect no Na D absorption in the Milky Way or 
LMC to a equivalent width limit of ~0.4 A. There is the possibility 
of Na D absorption in the host galaxy with a redshifted velocity of 
200 km s _1 relative to the host galaxy with an equivalent width of 
0.5 A. This possibility is represented by the red dashed lines. 



the host galaxy. There are three possible explanations for 
the lines. First, the lines could be a spurious detection 
caused by random noise. This is plausible given the level 
of noise in the spectrum and the slightly odd line profiles. 
Second, there could be an absorbing cloud traveling to- 
ward SN 2003ma. SN 2003ma is at the same redshift as 
its host galaxy; therefore, SN 2003ma (either in a disk or 
the halo) would have its velocity vector pointing mostly 
perpendicular to our line of sight. If the cloud is in the 
disk, then SN 2003ma must be in the halo, and likewise, 
if SN 2003ma is in the disk, then the cloud must be in 
the halo. Third, there is local material falling inward to- 
wards SN 2003ma. Either circumstellar material or gas 
from an H II region may be under gravitational collapse, 
however the velocity appears to be very high for such 
a scenario. Given the possibilities, we believe the most 
plausible is that there is no significant Na D absorption 
near the redshift of the host galaxy, which indicates no 
additional extinction from the host galaxy. However, the 
fact that there is a strong IR echo (see Section I4.8[) in 
combination with the type of the host galaxy being a 
star-forming galaxy which are generally dust-rich, makes 
it very likely that there is significant extinction in the 
host. The lack of Na D absorption might then be due 
to bleaching of the Na D doublet out to a few hundred 
pc by photoionizing Na I to Na II. Considering all these 
circumstances we assume as a lower limit no additional 
extinction from the host galaxy throughout the paper. 
Consequently, the total assumed E(B — V) = 0.348 mag 
is also a lower limit. Any additional extinction beyond 
the Milky Way Galaxy and LMC foreground component 
increases the required SN luminosity and hence makes 
the SN even more extreme in its characteristics. 

For each passband, we compute the combined Galactic 
and LMC extinction at its effective wavelength A e fi with 
A(Aeff) = {A(X cS )/A v } E(B - V)R V using Equation 1 



from iCardelli et al.l (|1989f ). Table [5] shows the effective 
wavelengths and extinction for the various passbands in 
the second and third column, respectively. 

4.2. Temperature at Peak 

As described in Scction f3.il we obtained a spectrum of 
the event and host galaxy on 2003 December 22, ~3 days 
after peak. To examine the temperature of this spec- 
trum, we focus on the dereddened, galaxy-subtracted 
spectrum (see Section 13.11 and Figure [7]) . We ignore the 
region of the spectrum corresponding to 6350 - 6800 A in 
our analysis below to avoid any Ha contamination. The 
spectrum shows a blue continuum without any signifi- 
cant emission or absorption lines. This is similar to the 
spectra of SNe II caught within days of explosion which 
also show a relative featureless blue continuum which 
can be fitted wit h blackbody temperatures of ~ 10,000 
- 20 000 K (e.g.. iDessart fc Hillier| [200l iDessart et all 
I2008t iGezari et all l2009al ). The earliest spectrum of 
SN 2005ip also demonstrates a blue, featureless contin- 
uum on day 1 after discovery that is fitted with a black- 
body with Tbb = 7300 K; however, this temperature is 
most likely a lower-limit due to the uncertainty in the 
amount of internal extinction ([Smith et aLll2009b|) . 

Figure [10] shows the 2003 galaxy-subtracted spectrum 
of SN 2003ma (black line), in comparison to black bod- 
ies of 5000, 10,000, 15,000, and 20,000 K. In order to 
further constrain the blackbody fit, we use the observed 
broadband fluxes in V7?sm and /ogle (blue square and 
red diamond). The horizontal error bars indicate the 
width of the filters. Black bodies with T < 10, 000 K 
do not fit well (see solid orange line and dashed sky 
blue). For higher temperatures, a differentiation is not 
possible since the observed wavelength range only covers 
the Rayleigh- Jeans tail of the blackbody. Therefore we 
can only set a lower limit of the temperature at peak of 
Tpcak > 15, 000 K. 

4.3. Absolute Magnitudes 

We calculate the absolute magnitude Mx in the rest- 
frame filter X from an observed magnitude my in filter 
Y as 

M x = my - A Y - DM - K X y (1) 

where Ay is the extinction in filter Y as described in 
Section 14.11 DM is the distance modulus and Kxy is 
the K-correction. For a redshift of z = 0.289, we com- 
pute DM = 40.86 mag. At this redshift, the filters B, 
VR, and I in the observers frame are best matched by 
the rest-frame filters U, B, and R, respectively. The 
K-correction Kxy as well as the bolometric correction 
described in the next section depend on the choice of the 
spectral energy distribution (SED) used to represent the 
SN SED. Due to the limited wavelength range of the ob- 
served spectra in addition to the imperfect background 
subtraction in particular for the later time spectra, we 
cannot use the SN 2003ma spectra to directly calculate 
the bolometric correction. Thus, we choose the following 
two SED models: 

• Solar model: The solar spectrum is often used 
to estimate the bolometric correction. Although it 
does not necessarily represent the true underlying 
SED, it is commonly used in the literature for other 
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Fig. 10.— Spectrum of SN 2003ma (black line) from 2003 Decem- 
ber 22 (~3 days after p eak) with galaxy contamination subtracted 
as discussed in Section 13.11 and dereddened. The galaxy Ho line 
does not fully subtract. Therefore the spectrum in the rest frame 
wavelength range 6350 - 6800 A is not reliable and we do not use 
it in the analysis. Overplotted are the dereddened observed fluxes 
derived from the VRsm (blue squares) and -/ogle ( re d diamonds) 
photometry. The horizontal error bars indicate the width of the 
filters. The solid orange line, dashed sky blue line, dot-dashed 
purple line, and dotted cyan line indicate blackbody fits of 5000, 
10,000, 15,000, and 20,000 K to the spectrum and photometry, 
respectively. 



SNe, and thus it allows for a direct comparison to 
them. 

• BB15K/88Z model: As we show in Section l4~2l 
the peak temperature is T pea k > 15, 000 K. Thus we 
choose a 15,000 K blackbody as our SED during the 
peak (MJD < 53010). Motivated by the fact that 
the light curves and color evolution of SN 2003ma 
are very similar to SN 1988Z after peak, we use 
the 6 April 1989 spectrum of SN 1988Z as the 
SED for the first plateau (53010 < MJD < 53110), 
and an average of three SN 1988Z spectra (cor- 
responding to dates 18 April 1990, 20 February 
1991, and 4 February 1992) for the second plateau 
(MJD > 53110). The particular SN 1988Z spec- 
tra were chosen for their large wavelength coverage 
(4500 - 9000 A) and high S/N. 

While the solar SED is simple and should provide reason- 
able measurements, the BB15K/88Z model is the more 
physical model since it attempts to accurately describe 
the SN SED. The fifth column in Tabic [5] and the third 
through fifth columns in Table [5] show the K-corrections 
for the two models. Note that the differences between 
the K-corrections derived from the different models are 
a few tenths of a magnitude. Figure [TT] and Q21 show 
the absolute magnitude light curve of SN 2003ma (blue 
squares and red circles) derived assuming the solar SED 
for the K-corrections. 

4.4. Luminosity & Energy 

We derive the luminosity Lx of SN 2003ma in a pass- 
band X and its bolomctric luminosity L by comparing it 



to the sun with 

L x =£ QiJC l0-°- 4 ( M *- M a.*) (2) 

Lho\,x = bxLx, (3) 

where bx is the bolomctric correction for band X. 
The solar luminosities L@ t x are c alculated with syn- 
thetic photometry off the solar SED (jLeieune et al.|[l997l) 
through Bessell passb ands, and th e solar magnitudes 
Mq } x are taken from iColina et al.l (|1996f ) . They are 
shown in the 6th and 7th column in Table respectively. 
We combine Ism and /ogle in all rest-frame measure- 
ments since they match well to the same rest-frame pass- 
band. We calculate the bolometric correction bx for a 
filter X as the ratio of the bolometric luminosity of the 
reference SED with its luminosity in X, where the ref- 
erence SED is defined by the model used (either solar 
model or BB15K/88Z model) and for the BB15K/88Z 
model also by its light-curve phase (see 8th column in 
Tableland 6th - 8th column in Table©. Note that 
the bolometric luminosity Lboi.x depends strongly on the 
bolomctric correction and thus on the filter X. Since we 
have good light curves in rest frame B and R, we can get 
estimates of the bolometric luminosity from both light 
curves and then compare them. In addition, we add the 
luminosities in the rest frame B and R together (Lb+r), 
and estimate the bolomctric luminosity Lboi.B+i?, using 
a bolometric correction bg , R = bg 1 + b^ 1 . Combining 
both filters decreases the systematic error since the over- 
all corrections are smaller, as we discuss in Section [4.51 
For both the solar model and the BB15K/88Z model 
we estimate the amount of energy emitted in the rest 
frame B, R, and the combined B + R passbands by 
integrating over time (see 10th and 9th column in Ta- 
bles [3] and HI respectively). The biggest source of sys- 
tematic error is caused by the gap in the data in the 
rest frame B (observed frame VR) between 53020 < 
MJD < 53250, which contains the transition between 
plateau 1 and plateau 2 (see Figure [4]). A simple lin- 
ear interpolation over the gap may introduce a system- 
atic overprcdiction of the flux. We estimate that this 
overestimate is smaller than 10%, and does not change 
our results. Note that the bolomctric correction using 
the BB15K/88Z model is only pseudo-bolometric, since 
the SN 1988Z spectra used only have a finite wavelength 
range from about 4000 - 10,000 A, not considering any 
contribution redward of the z band or blueward of the B 
band. 

4.5. Energy Emitted 

The upper left panel of Figure [14] shows the luminosi- 
ties in the rest frame B, R and combined B + R pass- 
bands using the solar SED for the determination of the 
K-correction. In the upper right panel, the same lu- 
minosities are shown, but using a 15,000 K blackbody 
SED for the K-correction during peak, and SEDs from 
SN 1988Z for the K-corrections during the plateau phase. 
The K-corrcctions for the two models are listed in Ta- 
ble [S] and [5] Note that the differences in K-corrections 
and thus luminosity are not bigger than 20%. 

We estimate the total emitted energy in the differ- 
ent passbands by integrating over time (see solid lines 
in lower panels of Figure [T4l and 10th and 9th column in 
Tables [5] and [H respectively). For the solar model, the 
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Rest-Frame Days Since Maximum 

Fig. 11. — Light curve of SN 20 03ma (DM = 40.86 m ag, z = 0.289) in absolute, rest-frame B and R band ma gnitudes in comparison to 
SN Iln 1988Z (DM = 34.77 mag ; lAretxaga et aI.llT999T), SN Hn 2005ip (D M = 32.36 mag; ISmithc t al. 2009J3), the e xtremely luminous 
SN Iln 2006gy (DM = 34.32 mag; ISmith et alJ20 07b; Agnolctto ct al. 200 <3), SN ILL 2008e s (DM = 40.08 mag, z = 0.21; IMiller et al.ll2009l ; 
IGezari et alJ2009a» . and SN IIL 2005ap ( DM = 40.81 mag, z = 0.283; lQuTmbv et al.l2007f l. SN la 1998aq (DM = 31.696,|Ricss ct al.j200l), 
and SN IIP 1999em (DM = 29.569 mag; ILeonard efal]|2002l) in rest-frame days since maximum. The times of maxima for SNe 1998aq, 
1999em, 2006gy, 2008es, and 2005ap are well constrained, but for SNe 1988Z and 2005ip there are no upper limits before the date of 
discovery, and so we assume that the date of discovery is the peak date. The light curves are in rest-frame B, R or i?-likc unfiltercd (UF) 
magnitude^. We do not apply any K-correction to the SNe with distances < 100 Mpc (SNe 1998aq, 1988Z, 1999em, 2005ip, and 2006gy). 
We convert the light curve of SN 2008cs from observed i band to rest-frame R band by applying a K-correction of —0.043 mag. We convert 
the light curve of SN 2005ap from unfiltercd magnitudes to rest-frame R band by applying a K-correction of 0.054 mag, assuming that the 
observed magnitudes are in R band. We apply extinction corrections of 0.036, 0.24, 0.126 and 1.78 mag to the light curves of SNe 1998aq, 
1999cm, 2006gy and 2005ip, respectively. The published light curves of SNe 1988Z, 2008es, and 2005ap arc already extinction corrected. 
The observed VR and /-band measurements of SN 2003ma are extinction corrected by 0.66 and 1.01 mag, respectively, and then converted 
into rest-frame magnitudes in B and R bands with K-corrections of —0.63 and —0.51 mag, respectively. All K-corrections are calculated 
assuming a solar spectrum. 

a Thc upper limits on the brightness of SN 2006gy are R band measurements from [Agnolctto ct al. (2009) 
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Fig. 12. — Light curve of SN 2003ma in absolute magnitudes in 
comparison to SN 1988Z in days since explosion, on a log scale. 
The time axis is in the rest-frame. Offsets have been added to the 
SN 1988Z absolute magnitudes in order to compare the shape of 
the light curves. Since the time of explosion is not well constrained 
for SN 1988Z, we shift the light curves by 12 days to match the 
bend in the I band light curve of SN 2003ma seen after 100 days. 



Fig. 13. — Rest-frame B — R colors of SN 2003ma in comparison 
to SN 1988Z in days since the estimated time of explosion to . The 
time of explosion of SN 1988Z is assumed to be 12 days before 
discovery. The time axis is in the rest-frame. Both SNe show an 
increase in B — R color until ~1000 days, after which the B — R 
color plateaus. 



emitted energy in rest-frame R band is 4.5 x 10 ergs 
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TABLE 5 



Observed 
Passband Y 


(A) 


A Y 
(mag) 


Restframe 
Passband X 


K X Y 
(mag) 


(10 32 ergs s _1 ) 


(mag) 


bx 


E x 
(10 50 ergs) 


Eho\,X 
(10 50 ergs) 


BsM 

VRsm 
Ism 

IoGLE 
VR + I 


4379 
5842 
7898 
7980 


1.44 
1.01 
0.66 
0.65 


U 
B 
R 
R 
B+R 


0.84 
-0.63 
-0.49 
-0.51 


1.89 
4.41 

6.48 


5.60 
5.47 

4.46 


20.40 
8.72 

5.93 

3.53 


4.76 (0.12) 
4.48 (0.10) 
9.24 (0.14) 


41.47 (1.06) 
26.57 (0.60) 
32.62 (0.49) 



Note. — Table of K-corrections and bolomctric corrections for the solar model. For each of the passbands Y (first column), 
the effective wavelength Xy , extinction Ay , the passband X that matches best in the restframe, the K-correction Kxyi the 
luminosity Lq ; x of the sun in X, the absolute magnitude Mq^x of the sun in X, and the bolometric correction bx are 
shown in the 2nd-8th column, respectively. The K-correction and bolometric correction are calculated using the solar SED. 
We combine Tsm and /ogle in all restframe measurements since they match well to the same restframe passband. The 
9th column shows the integrated energy Ex in passband X, and the 10th column the integrated bolometric energy -Ebol.x 
derived from the passband X photometry. The uncertainties quoted in Ex and i^bol X are the random uncertainties and 
don't contain any of the systematic errors due to e.g. the interpolation of the lightcurves, K-correction, etc. 



TABLE 6 



Observed 
Passband Y 


Restframe 
Passband X 


Peak 


Kxy 
Plateau 1 


Plateau 2 


Peak 


bx 
Plateau 1 


Plateau 2 


Ex 
(10 50 ergs) 


(10 50 ergs) 


VR SM 

Ism 
^ogle 
VR + I 


B 
R 
R 
B+R 


-0.44 
-0.57 
-0.57 


-0.54 
-0.59 
-0.62 


-0.61 
-0.74 
-0.80 


7.62 
12.22 
4.69 


26.66 
3.88 
3.38 


27.57 
6.15 
5.03 


5.04 (0.12) 
3.69 (0.08) 
8.67 (0.13) 


124.96 (3.40) 
22.28 (0.49) 
40.15 (0.63) 



Note. — Table of K-corrcctions and bolomctric corrections for the BB15K/88Z model. For each of the passbands Y (first column), 
the passband X that matches best in the restframe (second column), the K-corrcction KxYj an d the bolomctric correction bx arc 
shown. For both the K-correction and bolomctric correction, wc use 3 different SED's for the 3 different light curve phases peak, 
plateau 1, and plateau 2. For the peak, we use a 15,000 K blackbody SED. Since the lightcurve and color evolution of the first and 
second plateau is similar to that of SN 198SZ, we use a spectrum of SN 1988Z from 04/06/1989 as SED for plateau 1, and an average 
of SN 1988Z spectra from 04/18/1990, 02/20/91, and 02/04/92 for plateau 2. The K-correction for the peak, plateau 1 and plateau 2 
is shown in 3rd, 4th, and 5th column, respectively. The bolomctric correction for the peak, plateau 1 and plateau 2 is shown in the 
6th, 7th, and 8th column, respectively. We combine /sm and /ogle in all restframe measurements since they match well to the same 
restframe passband. The 9th column shows the integrated energy Ex in passband X, and the 10th column the integrated bolometric 
energy Ebol,X derived from the passband X photometry. The uncertainties quoted in Ex and E^ i f x are the random uncertainties 
and don't contain any of the systematic errors due to e.g. the interpolation of the lightcurves, K-corrcction, etc. 



with a similar amount emitted in the rest-frame B band 
of 4.8 x 10 50 ergs. This is more than ten times the en- 
ergy emitted by SN 1988Z in the B band over roughly 
the same period of ti me (~4.3 year) of ^3 x 10 49 ergs 
(jAretxaga et al.l[l999( ). By adding the luminosities in the 
rest frame B and R together, and integrating over the 
light curve, wc calculate the energy emitted in the com- 
bined rest frame B and R to be Eb+r = 9.2 x 10 50 ergs 
(black line in lower left panel of Figure [14)) . When we 
use the BB15K/88Z model for the K-correction, we get 
energies in B, R, and B + R of 5.0 x 10 50 , 3.7 x 10 50 , and 
8.7 x 10 50 ergs, respectively. These energies differ by only 
< 15% from the energies based on the solar model, and 
we thus estimate that our systematic errors are < 15%. 

The bolometric luminosities and hence bolomctric en- 
ergies have a strong dependency on the SED used to 
calculate the bolometric corrections. The dashed lines 
in the lower panels of Figure [14] show the bolometric 
energies derived from the different passbands using the 
solar and BB15K/88Z model on the left and right, respec- 
tively. The bolometric energies based on the solar model 
agree reasonably well and range between 2.7 x 10 51 and 
4.1 x 10 51 ergs. However, the bolometric energies derived 
using the BB15K/88Z model significantly differ, e.g., a 
factor of five between rest frame R (2.2 x 10 51 ergs) and 
B (12.5 x 10 51 ergs), even though this model is physically 
more motivated. 



This large difference between the rest frame B and 
R derived bolometric energy is mainly driven by the 
bolometric correction based on the SN 1988Z spectra 
for the two plateau phases. The late-time spectra of 
SN 1988Z-like SN are dominated by the strong Ha emis- 
sion line on the one hand, and the blue pseudo-continuum 
on the other hand. At late phases, this blue contin- 
uum is detected in most interacting SNe, e. g., SN 2006.jc 
( Foley et al.ll2007t ). and has been shown by iSmith et al.l 
( 2009bf ) to be fluorescence of a forest of emission lines 
that can only be resolved if the lines are narrow enough, 
e.g., SN 2005ip. The main difference between the spec- 
tral shape of these late-time spectra of SNe Iln is the 
relative strength of the Ha emission line (rest frame R) 
compared to the blue pseudo-continuum (rest frame B 
and bluer). The rest frame B - R color of SN 1988Z is 
significantly redder than the one of SN 2003ma, indicat- 
ing that for SN 2003ma the relative strength of the Ha 
line compared to the UV pseudo-continuum is weaker 
than for SN 1988Z. Thus using the spectra of SN 1988Z 
leads to a bolometric over-correction in the rest frame B 
band and an under-correction in the rest frame R band, 
which explains the large discrepancy between the two 
estimated bolometric luminosities. By combining rest 
frame B and R, this effect gets largely compensated, 
and the bolometric energy of 4.0 x 10 51 ergs is compa- 
rable to the bolometric energy using the solar model of 




Fig. 14. — Top left: The luminosities in the rest frame B, R and combined B + R passbands using the solar SED for the determination of 
the K-correction (solar model). Top Right: Same luminosities, but using a 15,000 K blackbody SED for the K-correction during peak, and 
SEDs from SN 1988Z for the K-corrections during the plateau phase (BB15K/88Z model). Bottom left and right: The integrated emitted 
energy in the B, R and combined B + R passbands using the solar an d BB 15K/88Z model, respectively, in solid lines. The dashed lines 
show the integrated bolometric energy for the two models. See Section 14.51 for a discussion. 



3.2 x 10 51 ergs. We thus estimate the bolometric radiated 
energy of SN 2003ma to be (3.6 ± 1) x 10 51 ergs, with the 
error dominated by the systematic error in the bolomet- 
ric correction. This is significantly more than the total 
radiated energ y of 2.3 - 2 5 xltP 1 ergs by SN 2006gy 
in -0.6 years (|Smith et al j |2009al). and 1.1 x 10 51 ergs 



for SN 2008es in ~0.4 years (jMiller et al.ll2009t) . and ten 
times the energy r adiated by SN 1988Z over 8.5 years 
of 3.2 x 10 50 ergs (jAretxaga et al.lll999T ). Analogous to 
SN 1988Z, the true bolometric luminosity of SN 2003ma 
is probably a factor of 10 higher if you include emission 
from the radio to X-rays that is produced from the in- 
teraction of the SN shock with the CSM, resulting in a 
total radiated energy budget of at least 10 52 ergs. 

Even though SN 2003ma is not quite as bright at peak 
as the most luminous SNe like SNc 2005ap, 2006gy, and 
2008es, its total emitted energy is larger since half of 
the radiated energy is emitted after 350 days due to its 
very long decline (see the dotted lines in Figure [14] for 
Eboi,B+n), at times when the other extremely luminous 
SNe do not produce any significant radiation. Such a 
large amount of radiated energy can be achieved if two 



things coincide: (1) the initial explosion has a large total 
(combined radiative and kinetic) energy and (2) most of 
the kinetic energy is converted into radiation by interac- 
tion with a CSM. 

4.6. Peak Luminosity 

The peak luminosity of a SN II from thermal energy 
deposited into the ejecta following shock break out can 
be estimated to be 



U = 



(3c E SN R 
2k M 



(4) 



= 2.5 x 10 4 



E 



1 ergs 



R 



10 



14, 



(5) 

where £?sn is the energy of the SN explosion, M is the 
mass of the ejecta, Rq is the initial radius of the star, 
(3 is a numerical constant related to the radiat ive diffu- 
sion time, and n is the opacity (jArnettl [l996j). Thus, 
the most luminous core-collapse SNe will result from 
progenitors with large initial radii and low densities. 
Given the measured total energy emitted in the optical 
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of -Bsn = 4 x 10 51 ergs, and L = 5 x 10 43 ergs s~\ this 
implies Ro = 2 x 10 14 \Mj (10M©)] cm, a value approach- 
ing th e radii of the largest known RSGs (jLevesque et al.l 
120051 ) . If the progenitor were an Eta Carinae-type mas- 
sive star with M ~ 100M©, this would require a radius 
ten times larger, which is difficult to reconcile with stellar 
radii predicted by stellar evolutionary models. 

If the peak of the light curve is powered solely by ra- 
dioactive decay, we expect 



M m /M Q 



L 



1.42 x 10 43 ergs s" 



a (t/llldays) 



(6) 



For SN 2003ma, which had i poak = 5 x 10 43 ergs s" 1 at 
t = 8.5 days after explosion in the SN rest-frame, the ini- 
tial 56 Ni mass is 4 M Q , which is one to two orders of mag- 
nitude larger than typically produced in SNe IIP (jHamuvl 
l2003h . Although this amount of 56 Ni could be produce d 
in a pair-instability explosion (|Scannapieco et aLll2005l ). 
the slow rise to maximum and constant slope of decay af- 
ter maximum in such models are incompatible with the 
light curve of SN 2003ma. 

4.7. Kinetic Energy 

Measuring the kinetic energy from a SN typically re- 
quires an cjecta mass and a velocity for the ejecta. For 
SN 2003ma we were unable to measure the velocity of the 
ejecta from any of our spectra. However, we can estimate 
the energy converted from shock energy to kinetic energy 
by examining the velocity of the post-shock material. If 
we assume that the pre-shock circumstellar material has 
a coherent velocity, v w and a mass, m, it has a kinetic 
energy, 



(7) 



If we assume that the post-shock material has a one- 
dimensional RMS velocity, w ps , then the post-shock ki- 
netic energy of the circumstellar material is 



and the shock deposited an energy, 

AE = im(3t& - vl). 



(8) 



(9) 



If we assume that the pre-shock circumstellar material 
has a constant density, /?, and that the shock traveled at 
a velocity of v s with no appreciable deceleration, then 



m = -itp{v s t) 



(10) 



is the mass of circumstellar material swept up by the 
shock in a time t after explosion. 

Combining Equations [9] and [10] and assuming that 



w ps > w w , we find 



AE = 2npvl a v 3 s t 3 . (11) 

If we have a constant conversion of kinetic energy to lumi- 
nosity, then L = r]AE/t, where r\ is the efficiency factor. 



The density can then be written as 
1 L 



p = 5.6x 10~ 20 - 



rj 10 42 ergs s" 1 V6000 km s" 

-3 / ± \ -2 



■ ps 



20,000 km s~ 



t 



year 



g cm 



(12) 



^From Equations [TU] and Q2] we have the total swept-up 
mass, 



0.03- 



L 



n 10 42 ergs s- 1 V6000 km s 



-2 



t 



year 



M e , 
(13) 

which is independent of any assumed ejecta velocity. At 
the time of the last spectrum, when L = 5 x 10 42 ergs s _1 , 
v ps = 6000 kms" 1 , and t = 3.9 years, SN 2003ma had 
swept up 0.6 Mq of circumstellar material and had at 
least 6.4 x 10 50 ergs of kinetic energy if 77 = 1. These 
values are lower limits since if r\ < 1, the swept-up mass 
and kinetic energy increases. 

4.8. Dust 

The Spitzer photometry in IRAC bands 1-4 and MIPS 
band 1 (see Table \2§ is shown in Figure [TS] with respect 
to the rest wavelengths (2.8, 3.5, 4.5, 6.2, and 18.5 pm, 
respectively) corresponding to the observed wavelengths 
of 3.6, 4.5, 5.8, 8.0, and 24 pm, respectively. There is a 
significant difference between the first epoch (2005 July 
15 - 26; ^580 days after maximum) and the second epoch 
(2005 October 26 - November 2; ^675 days after maxi- 
mum) at restframe 3.5, 4.5, and 6.2 /im: the fluxes at 3.5 
and 6.2 pm increase by 4er and 10cr, respectively, whereas 
the 4.5 pm flux decreases by 4<r. Note that we use the 
relative flux error for the comparison of fluxes at differ- 
ent epochs, but the same filter (4th column in Tabled]), 
as we have discussed in Section [2.21 For all comparisons 
of fluxes from different bands, we use the total errors 
(6th column in Table [2]) , which includes the wavelength 
dependent error introduced by the crowdedness of the 
field. 

Such a significant change in flux over such a short 
time span makes it highly likely that at least some of 
the measured infrared flux is caused by the SN. Note, 
however, that the direction of these differences across 
the bands is not in the same direction in neighboring 
bands, and therefore SN 2003ma would not have been 
considered a variable candid ate according to the crite- 
ria used by iViih et al.l pOOl ) for the SAGE-LMC vari- 
ables list. The observed 24 pm flux is essentially un- 
changed and is either an IR echo as observed for other 
SNe (e.g., SN 2 002hh. IBarlow et all [2001 IMeikle etall 
12001 SN 2004et. lKotak et aLteOOgfi or from other sources 
in the host galaxy. 

The changes in flux cannot be explained by a sim- 
ple change in temperature and/or flux of a black- 
body, but rather requires emission bands of molecules 
or dust. The first detection of carbon monoxide 
(CO) in a SN in the fundamental (4.6 pm) and 
first overt one (2.3 pm) bands was in the spectra of 
SN 1987A (ICatchpole et"aIIU987t IMcGreeor et alJ ll987t 
lOliva et al.l 119871: iSpyromilio et al.l 11988ft . In~~addition 
to CO, silicon monoxide (SiO) emission ba nds were de- 
tected in SN 1987A at 7.9 and 8.5 pm ([Roche et all 
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U99l IMeikle et all U993 iLiu fc Dalgarnol U99l . Sub- 
sequently, MIR CO lines were d etected in SNe 1995ad 
(ISpyromilio fc Leibundg ut] 1199611 1998S dGerardy et al.l 
l2000 HFassia et al.ll200lD."l998dl (ISp yromilio et alJl200lH 
1999em (IS pyromilio et aLll2001ft. 2000ew dGerardv et al.1 
I2002H ) and 2004et (|Kotak et all 12009( 1 . which also 
showed SiO lines. The wavelengths of the CO and SiO 
emission bands are indicated in Figure 1151 The match 
between the CO and SiO lines with the observed flux 
changes is poor, except for possible contribution of the 
CO fundamental band at 4.6 /im. However, we note that 
for both SNe 1987A and 2004et the CO and SiO bands 
appeared ~100 days after maximum and disappeared by 
~500 days, significantly earlier than the increase seen in 
the second epoch for SN 2003ma. 

Another possible explanation is that the strong and 
variable emission detected in the 8 /jm IRAC band is 
from the 6.2 /xm emission feature of polycyclic aromatic 
hydrocarbon (PAH) nanoparticles. Such emission is 
routinely detected in the Mi l ky W a y and other galax- 
ies (e.g., IGenzel fc Cesarskvl l2000l; IVoeler et all 120051: 



Irwin fc Madden! 120061 : iDraine fe Lil 120071 : iSmith et al.l 
2007all . 
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Fig. 15.— Spitzer IRAC (3.6, 4.5, 6.0, and 8.0 (im) and MIPS 
(24 fim) photometry versus rest wavelength, for 584 and 683 days 
after maximum. The solid curves show models which include 
hot dust, warm dust (see text), and a stellar component with 
a Raylcigh-Jeans spectrum. Filled diamonds are the model con- 
volved with the IRAC and MIPs response functions. The photom- 
etry is consistent with the observed emission consisting of steady 
emission from warm dust in the galaxy, plus a time-variable iR 
echo produced by hot dust. 



Here we present a dust mo 
observed IR emission. First, 
of the emission observed at 3 
is stellar continuum and we 
wavelengths assuming a Raylc 
dashed line in Figure [T5|) . To 
a mixture of PAHs, graphite, 



del that could explain the 
we assume that the bulk 
.6 /im (rest frame 2.3 /im) 
extrapolate this to longer 
igh- Jeans spectrum (short- 
this we add emission from 
and silicate dust, from the 



dust models of IDraine fcTH (|2007f l. For SN 2003ma we 
assume the radiation field heating the dust to have the 
spectrum of a 15,000 K blackbody, and have solved for 
the temperature distribution functions for different grain 
sizes and compositions, and different radiation intensi- 
ties. The emission model has two components: steady 
emission from the "warm" dust in the galaxy, plus a 
time-varying IR echo from "hot" dust heated by radi- 
ation from the SN. 

To reproduce the 24 /im photometry, the warm dust 
must be hot enough for the silicate 18 /im feature to 
emit strongly. This is accomplished by a radiation field 
with U « 10 5 5 (U is the dust heating rate relative to 
heating by the starlight background in the solar neigh- 
borhood). This dust component is assumed to have the 
same ionization fraction as a function of grain size as as- 
sumed by DL07 for the Galaxy, and which has been found 
to work well for the IR-subm m emission from, e.g ., the 
galaxies in the SINGS sample (jDraine et all 12007(1 . The 
PAH abundance in the warm dust component must be 
low, gpAH ~ 0.5%, in order to not overproduce the ob- 
served flux in IRAC bands 2 and 4. This component has 
a dust mass of ~25OM ; for a normal dust-to-gas ratio, 
this corresponds to a hydrogen mass of ^2.5 x 10 4 M Q . 
This is a plausible amount of gas for hot photodissocia- 
tion regions in a starburst galaxy. The spectrum of this 
emission is shown in Figure 1151 The total IR luminosity 
of this dust component is ~1.1 x lO lo i . There is pre- 
sumably additional emission from a much larger mass of 
cool dust in the galaxy, but we lack longer wavelength 
observations to constrain this component. 

The dust producing the IR echo is assumed to be very 
strongly heated, with U « 10 7 5 . The IR echo dust mass 
is ~0.3M Q on day 584, and ~1A/ Q on day 683. The 
total IR luminosity of this component on day 683 is ~ 
4 x 1O 9 L , requiring that the IR echo dust be absorbing a 
substantial fraction of the SN light. This dust is assumed 
to have a PAH abundance <7pah = 4.7%, similar to the 
PAH abundance in many solar-mctallicity star-forming 
galaxies. The PAHs are assumed to be ionized by the 
radiation from the SN; the PAH neutral fraction in this 
component is assumed to be negligible (the absence of 
neutral PAHs lowers the ratio of 3.3 /jm emission relative 
to 7.7 and 8.6 emission, which is needed in order to 
not overproduce the emission into IRAC band 2). The 
intense radiation heats the large grains to T sa 300 K; 
the contribution of the broad 9.7 /im silicate feature is 
apparent in this component of the model spectrum. Even 
in this intense radiation field, the emission in the 3.3 and 
6.2 /im PAH features continues to be due primarily to 
single-photon heating. 

This two component model, with the mass of the hot 
dust increasing by a factor ~3 between days 584 and 
683, reproduces most of the IR photometry. There is 
admittedly a discrepancy with the observed flux density 
in IRAC band 3 (A rost = 4.5 /im): the observed flux de- 
creases by a factor ^3 between days 584 and 683, and our 
model does not reproduce this. Perhaps CO 1-0 4.6 /tm 
emission (not included in our model) contributes to the 
high flux measured on day 584. 

The dust heating parameter U = 10 7 5 corresponds to 
an energy density u(hv < 13.6eV) — Ax 10 -5 erg cm -3 . 
The distance of the dust from the SN is not known, but 
it is reasonable to assume a distance ^600 light-days 
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(0.50 pc). For this distance, the SN luminosity required 
to produce U = 10 7 5 is ~1.4 x 10 10 L Q - consistent with 
the photometry in Figure [Ml 

A sphere of radius 0.5 pc with density nn = 8 x 
10 3 enr 3 contains Md w IMq of dust (assuming 
M d /M n ~ 0.01). Densities n H > 10 4 cm" 3 would be 
expected for cool gas in the high-pressure environment 
of a starburst galaxy. The IR echo requires ~ 1-M© of 
dust at a distance ~0.5 pc from the SN to produce the 
IR echo. This dust could be located in an interstellar 
cloud, or possibly in a swept-up shell of interstellar mat- 
ter surrounding a stellar-wind bubble. 

To summarize: (1) PAH emission features can explain 
the increase in the observed 4.5 and 8 /im flux densi- 
ties; (2) the required dust mass of « IMq should be 
~0.5 pc from the SN, which is reasonable; (3) the inten- 
sity of the SN radiation at this distance is consistent with 
what is required to heat the dust to reproduce the ob- 
served emission. We conclude that the observed IR pho- 
tometry of LMC 902 is consistent with an IR echo pro- 
duced by normal interstellar dust at a distance ~0.5 pc 
from the SN. We note, however, that the data does not 
sufficiently constrain the model to make it unique, and 
other explanations are possible. 



5. CONCLUSIONS 

This paper presents photometry and spectroscopy of 
SN 2003ma, an exceptional SN Iln at a redshift of z = 
0.289 behind the Large Magellanic Cloud. 

SN 2003ma has an /-band light curve with two long 
plateaus, and strong emission in the V and I bands for 
^1000 days, reminiscent of the long-lived emission seen 
in interacting SNe Iln like the prototype of the class, 
SN 1988Z. However, in contrast to all other SNe Iln with 
long-lived plateaus, SN 2003ma also has an exceptionally 
bright peak luminosity of = —21.5 mag, placing it in 
the class of extremely luminous core-collapse SNe, such 
as SNe 2005ap, 2006gy, 2008es, and 2008fz. The combi- 
nation of high peak luminosity and long-lasting plateau 
makes SN 2003ma more luminous than any other SN at 
80 days past maximum brightness. Consequently, half of 
the radiated energy is emitted after 350 days past maxi- 
mum light, while other extremely luminous SNe produce 
no significant radiation at such late times. 

The derived bolometric luminosity, and hence bolomet- 
ric energy, has a strong dependency on the SED used to 
calculate the bolometric corrections. In order to esti- 
mate the systematic bias of the assumed SED, we use 
different SEDs to calculate the bolometric corrections. 
Integrating the bolometric light curve, we find that a to- 
tal bolometric output of ^4 x 10 51 ergs was emitted over 
4.7 years - more than twice the total radiated energy of 
SNe 2006gy or 2008es, and ten times the energy radiated 
by SN 1988Z, making SN 2003ma the most energetic SN 
ever observed with respect to its radiative output. This 
value alone is within an order of magnitude of the max- 
imum energy that a core-collapse explosion can produce 
outside of the neutrino channel. 



The nearly flat late-time light curve over the long time- 
span of several years in combination with the persistent 
single-peaked intermediate-width Ha emission line are 
best explained by the conversion of the kinetic energy of 
the explosion into radiation by strong interaction of the 
SN with a dense circumstellar material. SN 2003ma is 
distinct from proto-typical SNe Iln in that the interme- 
diate width line is unusually broad, and it does not have 
a broad component. If SN 2003ma has narrow emission 
lines cannot be answered conclusively since it is difficult 
to separate the narrow emission lines from the SN, if they 
exist, from the narrow emission lines of the starbursting 
host. 

The fast velocity measured for the intermediate- width 
Ha component (^6000 kms -1 ) also points towards an 
extremely energetic explosion (> 10 52 ergs) which im- 
parts a faster blast-wave speed to the postshock mate- 
rial and a higher luminosity from the interaction than 
is observed in typical SNe Iln such as SN 1988Z. The 
large amount of radiated energy can only be achieved if 
the initial explosion has a large total radiative and ki- 
netic energy and most of its kinetic energy is converted 
into radiation by interaction with a CSM. A conventional 
core-collapse SN explosion can produce at most a few 
times 10 52 ergs of energy that can couple to baryonic 
matter. The measured kinetic energy of SN 2003ma is 
comparable to this limit (and similar to GRB-associated 
SNe) , and therefore SN 2003ma may have been produced 
by an alternative explosion mechanism such as the pair- 
instability SN. 

The Spitzer IR photometry shows significant variations 
in the IR from 584 days to 683 days after maximum in 
rcstframe 3.5, 4.5, and 6.2 /urn. Such a significant change 
in flux over such a short time span makes it highly likely 
that at least some of the measured infrared flux is caused 
by the SN. The changes in flux cannot be explained by a 
simple change in temperature and/or flux of a blackbody, 
but rather requires emission bands of molecules or dust. 
We find that the observed increase in the IR flux density 
at 4.5 and 8.0 /im is consistent with an IR echo of the 
SN produced by normal interstellar dust at a distance 
~0.5 pc from the SN in combination with PAH emission 
features. 

By any measure, SN 2003ma is a remarkable ob- 
ject. Its discovery and extensive lightcurve coverage 
were made possible by long-duration microlcnsing sur- 
veys. Its prompt identification and initial spectrocopy 
reveal the power of well-implemented reduction pipelines 
and planned, timely coordination of large telescope fa- 
cilities. The inclusion of satellite data provided im- 
portant constraints of the environment of this object. 
Ground-based facilities devoted entirely to deep, time- 
domain discovery such as the Large Synoptic Survey 
Telescope and PanSTARRS are almost certain to re- 
veal many more such objects and lead us towards even 
greater understanding of the progenitors required to pro- 
duce extremely luminous core-collapse supernovac like 
SN 2003ma and illuminate the full range of observational 
characteristics for what, for now, remains a small and 
elusive set of truly cataclysmic events. 
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APPENDIX 
CLASSIFICATION OF SN 2003mA 
In this section we summarize why we classify SN 2003ma as a SN and not one of these other options. 

Event in Milky Way or LMC 

From Figure^ we sec that SN 2003ma is 0.064±0.012 arcsec from a resolved source in our template images. In 2008, 
SN 2003ma was much fainter than this source; therefore, our 2008 spectrum is dominated by this source, a starburst 
galaxy at z = 0.289. If the event is in the LMC or Milky Way, then it must be either a transient or variable source 
with a quiescent brightness > 5 mag fainter than at peak. This excludes most flaring variable stars such as M dwarfs 
and Ac/Be stars since the amplitude of their flares are generally small compared to the total flux. Additionally, we 
would not expect a prolonged, > 5 year decline from peak. 

Another candidate is a nova, which could flare by > 5 mag. However, our spectra from 2003 does not show any 
emission lines at or near zero redshift, which we would expect since the event is 0.5 magnitude brighter than the 
source at that time and therefore should dominate the spectrum. Furthermore, the light curve is unlike any nova yet 
observed. 

Considering all this, it is unlikely that SN 2003ma occurred in either the LMC or Milky Way. 

Supermassive Black Hole Event 
AGN are often variable at various time-scales. Here we discuss if it is likely that SN 2003ma is an AGN. 



If SN 2003ma originates from a black hole in the supposed host galaxy then the supermassive black hole is offset 
by 500 pc from the center of the galaxy. 

We calculate the weighted average of 42 SuperMACHO VR difference image fluxes before the event (MJD 
< 52965) spanning 2 years, and find that it is within ler of zero, with a reduced x 2 = 1.2. None of the 42 
detections deviates from zero by more than 3cr, and typical 3<r lower limits on variability are 24.2 magnitudes. 
If we split up these 42 detections into three blocks 52224 < MJD < 52291, 52547 < MJD < 52643, and 
52932 < MJD < 52965, we find again that the average fluxes are consistent with zero with reduced x 2 of 1.01, 
0.96, and 1.73, respectively. Therefore the pre-event variability in VR is fainter than 24.2 magnitude during this 
2-year time period. 

We perform a similar analysis with the / band light curve and find that the pre-event OGLE-II (MJD < 52091) 
and OGLE-III (52091 < MJD < 52965) data is consistent with being constant with a reduced x 2 of 1.65 and 
1.51, respectively. We 3cr-clipped 4% of the data. When we bin the data (see Section [2~T1 and bottom panel in 
Figure [3]), we find that all average fluxes are within 3tr of zero, with reduced x 2 ranging from 0.67 to 1.91. 

Another test for the presence of AGN activity in the host galaxy are the ratios of the narrow emission lines, 



which are measured for all spectra an d presented in Table 
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in the standard diagnostic line ratios ((Baldwin et al.|[T981 " iVeilleux fc dsterbrocld[l987f ) classify the host galaxy 



as a star-forming galaxy. Although the emission lines may have a contribution from the transient object, in 
the post-peak spectra from 2004 — 2008 the [O III]/H/3, [N II]/Ha, and [S II] /Ha ratios are constant within the 
errors, which implies that at late times the narrow-line emission is dominated by the host galaxy. 

• We consider the rare case in which an otherwise dormant supermassive black hole tidally disrupts a star, and 
emits a flare of radiation from the accretion of the stellar debris. This could explain that the event itself is 
about 0.5 mag nitude brighter than the host and sh ows no var iability pre-eyent. These events a re typically bright 
in the X-rays (jKomossal l200l lEsauei et all l2007f ) and UV (jGezari et al.ll2006l 120081 l2009bl). however s everal 
candidates have emerged that have also been detected at optical wavelengths (jGezari et al.l 120081 . 12009b0 . The 
strong evolution of the foreground extinction-corrected colors of SN 2003ma from V — / « at peak to V — I > 1 
for times later than 1 year after the peak (see Section l4~3j) . however, are inconsistent with the steady blue 
colors expected for the blackbody sp ectral energy distrib utions and effective temperatures associated with tidal 
disruption events (T BB > 1 x 10 4 K; lGezari et"aLll2009bl ). 

The narrow-line ratios, listed in Table IH are all consistent with a star-forming galaxy. This, in combination with 
the lack of variability in the densely sampled SuperMACHO and OGLE difference imaging data over a baseline of 7 
years before the SN, place strong constraints on the presence of an AGN in the host galaxy. 

Supernova 

Since SN 2003ma is unlikely to be caused by a supermassive black hole or Galactic or LMC variable, the remaining 
likely transient phenomenon is that of a SN. Its blue continuum in 2003 (see Figures 151 and fT0 | and ^6000 kms" 1 
FWHM emission lines consistent with Ha in 2004 and 2008 are reminiscent of a SN of type Iln (lSchlegell H990). This 
identification is also consistent with the host galaxy being a star- forming galaxy. In Sections 13.11 and 13.21 we compare 
the spectroscopy and photometry of SN 2003ma to other SNe from the literature, and we find that the most consistent 
explanation for this event is that it is a SN Iln. 
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